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Abstract 

This paper analytically characterizes optimal sensor placements for target localization and tracking in 2D and 3D. Three types 
of sensors are considered: bearing-only, range-only, and received-signal-strength. The optimal placement problems of the three 
sensor types are formulated as an identical parameter optimization problem and consequently analyzed in a unified framework. 
Recently developed frame theory is applied to the optimality analysis. We prove necessary and sufficient conditions for optimal 
placements in 2D and 3D. A number of important analytical properties of optimal placements are further explored. In order to 
verify the analytical analysis, we present a gradient control law that can numerically construct generic optimal placements. 

Index Terms 

Fisher information matrix; gradient control; optimal sensor placement; target tracking; tight frame. 

I. Introduction 

Target localization and tracking using networked mobile sensor platforms has become an active research area in recent 
years. When localizing a target from noisy measurements of multiple sensors, the sensor placement can significantly affect the 
estimation accuracy of any localization algorithms. The term sensor placement as used here refers to the relative sensor-target 
geometry. This paper will address the optimal sensor placement that can minimize the target localization uncertainty. 

In the literature, there are generally two kinds of formulations for optimal sensor placement problems. One is optimal control 
|[l]-||4] and the other is parameter optimization ||5]-([T3|. 

The optimal control formulation is usually adopted for cooperative path planning |[T]-|[3|, where the aim is to estimate the 
target position on one hand and plan the path of sensor platforms to minimize the estimation uncertainty on the other hand. This 
problem is also referred as simultaneous-localization-and-planning (SLAP) |j2j. In a SLAP problem, target motion and sensor 
measurement models are considered as process and measurement models, respectively. The Kalman filter usually is applied 
to estimate the target position and to characterize the estimation covariance. In order to minimize the estimation covariance, 
an optimal control problem will be formulated. The disadvantage of this kind of formulation is that the optimal control with 
various constraints generally can only be solved by numerical methods. Analytical properties usually cannot be obtained. 

To avoid complicated optimal control problems, many studies including our work in this paper formulate the optimal 
placement problem as a parameter optimization problem. This kind of formulation has a long history and has been investigated 
extensively in p|-||T3|, to name a few. The parameter optimization formulation is based on the assumption that a rough 
estimation of the target position has akeady been obtained in other ways. By using this rough estimation, the sensor positions 
are the parameters to be optimized such that the consequent target localization based on the optimized placement will be more 
accurate. The objective function in the parameter optimization formulation usually involves the Fisher information matrix (FIM). 
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The FIM is the inverse of the Cramer Rao lower bound (CRLB), which is the minimum achievable estimation covariance. An 
unbiased estimator that achieves the CRLB is called efficient. An optimal placement, which maximizes a function (such as the 
determinant) of the FIM, can be interpreted as maximizing the target information gathered by the sensors or minimizing the 
estimation covariance of any efficient estimators. 

In contrast to the optimal control formulation, the parameter optimization formulation can be solved analytically. The 
analytical solutions are important for us to get insight to the effect of sensor placement on target localization uncertainty. It is 
notable that the numerical results based on the optimal control formulation generally are consistent with the analytical results 
based on the parameter optimization formulation. For example, the numerical simulations in |[T], ||2j show that the final optimal 
angle subtended by two sensors at the target is 90 degrees. In |j5), (8), Q, it is also analytically proved that the placement of 
two sensors is optimal if the angle subtended by the two sensors at the target is 90 degrees. 

In this paper, we will investigate optimal sensor placement by adopting the parameter optimization formulation. Our aim is to 
determine the optimal sensor-target geometry based on an initial estimation of the target position. Optimal sensor placement is 
of not only theoretical interest but also significantly practical value. Many studies have shown that target tracking performance 
can be improved when sensors are steered to form an optimal placement. In this paper, we only focus on determining optimal 
placements and will not address target tracking. One may refer to [9 1 for a comprehensive example that illustrates the application 
optimal sensor placements to cooperative target tracking. 

Until now, most of the existing studies only consider 2D optimal sensor placements p)-pO), p2) . Very few works in the 
literature have tackled 3D cases | [T3) . Analytical characterization of generic optimal sensor placements in 3D is still an open 
problem. In this paper we will extend the results in Q, ||7|-||9) from 2D to 3D. The extension will be non-trivial. Maximizing 
the determinant of the FIM has been widely adopted as the criterion for optimal placements in 2D. This criterion, however, 
cannot be directly applied to 3D cases because the determinant of the FIM is hardly analytically tractable in 3D. Motivated by 
this, we will propose a new criterion for optimal placement, which enables us to analytically characterize optimal placements 
in 2D and 3D. The existing analysis of 2D cases can be regarded as a special case of our general analysis for both 2D and 
3D cases. 

The existing work on optimal sensor placement has addressed many sensor types including bearing-only ||5), (8], 1 14 1, range- 
only (5), (9), IB), received-signal-strength (RSS) |j7|, time-of-arrival (TOA) (5), ||6|, and time-difference-of-aiTival (TDOA) 
Q, p2) . These sensor types are analyzed individually in the literature. One of the contributions of this paper is to unify the 
analyses of bearing-only, range-only, and RSS sensors. Based on our proposed optimality criterion, we will show the objective 
functions for the three sensor types are exactly the same. Hence their optimal placement can be analyzed in a unified way. 
Since the measurement models and FIMs of TOA and TDOA sensors are significantly different from those of the three, we 
will not consider TOA and TDOA sensors in this paper. 

By employing recently developed frame theory, we prove necessary and sufficient conditions for optimal placements in 2D 
and 3D. Frames provide a redundant and robust way for representing signals and are widely used in signal processing. We 
refer to p6) , p7) for an introduction to frames. It might be interesting to ask why frames arise in optimal sensor placement 
problems. This question can be loosely answered from the redundancy point of view. As mentioned in p6] , one would use 
frames when redundancy is a must. In our work, the redundancy can be interpreted as the ratio between the number of sensors 
and the space dimension. When the sensor number equals the dimension, there is no redundancy in the system, then we will 
show that the necessary and sufficient condition of optimal placement can be proved without using frames. But when the sensor 
number is larger than the space dimension, our optimality analysis will heavily rely on frame theory. 
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The paper is organized as follows. Section III] introduces preliminaries to frame theory. Section [III] presents a unified 



formulation for optimal placement problems of bearing-only, range-only, and RSS sensors in 2D and 3D. In Section IV 



necessary and sufficient conditions for optimal placement in 2D and 3D are proved. Section |V] further explores a number 



of important properties of optimal placements. In Section VI a gradient control law is proposed to numerically verify our 
analytical analysis. Conclusions are drawn in Section |VII| 



II. Preliminaries to Frame Theory 

Frames can be defined in any Hilbert space. Here we are only interested in d-dimensional Euclidean space M.'^. Let || • || be 
the Euclidean norm of a vector or the Frobenius norm of a matrix. As shown in p6)-p9), a set of vectors in M'* 

(n > d) is called a frame if there exist constants < a < 6 < +oo so that for all x G M'' 

n 

a\\x\\' <J2{x,^.f<b\\x\\\ (1) 

i=l 

where (•,•) denotes the inner product of two vectors. The constants a and b are called the frame bounds. A frame {(pi}'^^i 
is called unit-norm if = 1 for all i E {l,...,n}. Denote <i> — [(pi, . . . , (pn] G M''^". Because {x,(pi)'^ — [x^LpiY = 
x'^LPiipfx, inequality ([T]i can be rewritten as 

a\\xf < a;^$$^x < b\\xf, 

where the matrix = J27=i fiff called the frame operator. The frame bounds a and h obviously are the smallest and 
largest eigenvalues of ^i"^, respectively. Since a > 0, is positive definite. It is well known that d vectors in form a 
basis if the vectors span W^. Frame essentially is a generalization of the concept of basis. The frame {'Pi]'^^i can also span 
because "I>"I>^ is positive definite and hence $ is of full row rank. But compared to a basis, a frame have n — d redundant 
vectors. The constant n/d is referred as the redundancy of the system. When n/d= 1, the frame would degenerate to a basis 

of wK 

Tight frame is a very important concept in frame theory. A frame is called tight when a = b. From ([T]i, it is easy to see the 
frame {(pi}2^i is tight when 

n 

^fiVl^ald,. (2) 

1=1 

Taking trace on both sides of Q yields a = X]"=i ll'^'ilP/'^- is an important and fundamental problem in frame theory to 
construct a tight frame {(pi}f^i that solves (|2]) with specified norms. This problem is also recognized as notoriously difficult 
pO| . One approach to this problem is to characterize tight frames as the minimizers of the frame potential 

n n 

FP({^Jr=i)=EE(^^^^)'- ^3) 

i=l 3 = 1 

Frame potential is first proposed in |18j for unit-norm frames, and then generalized in | |T9) for frames with arbitrary norms. 

The following concept irregularity p6) , p9j is crucial for characterizing the minimizers of the frame potential. 

Definition 2.1 (Irregularity): For any positive non-increasing sequence {c;}"^]^ with ci > • • • > c„ > 0, and any integer d 
satisfying 1 < d < n, denote as the smallest nonnegative integer k for which 

1 " 

^ti < ^ E (4) 

i=k+l 
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The integer fco is called the irregularity of {c^}"^]^ with respect to d. 

Remark 2.2: The irregularity of a sequence is evaluated with respect to a particular positive integer The irregularity of a 
sequence may be different when evaluated with respect to different positive integers. In this paper, we will omit mentioning 
this integer when the context is clear. 

Because the index k — d — 1 always makes Q hold, the irregularity fcp always exists and satisfies 

< fco < c'- 1- 

When fco = 0, inequality (Hh degenerates to the fundamental inequality p9) 



1 

max < - y^c? (5) 
i=i,...,n ^ - d^ 

i—l 

In this paper we call the sequence {ci}f^^ regular when fco — 0, and irregular when 1 < ko < d — 1. The fundamental 
inequality (|5]l intuitively means that: a sequence is regular when no element is much larger than the others. Next we show 
several examples to illustrate the concept of irregularity. 

Example 2.3: Consider a sequence {ci}"^^ with c = ci = • • • = c„ and any d < n. Because l/dY^"^^ cf — nc^ /d > c^, 
the fundamental inequality (|5]l holds. Thus this sequence is regular with respect to any integer d < n. 

Example 2.4: Consider a sequence {ci}f^i = {10, 1, 1, 1} and d = 3. Note the feature of this sequence is that one element 
is much larger than the others. Because 10^ > 1/3(10^ + 1 + 1 + 1), the sequence is irregular with respect to d = 3. In order 
to determine the irregularity fco, we need to further check if {ci}f^2 = 1' 1} is regular with respect to d — 1 = 2. Since the 



elements of {ci}f^2 equal to each other, {ci}f^2 is regular with respect to 2 as shown in Example 2.3 Hence the irregularity 
of {ci}j^i with respect to d = 3 is fco = 1. This example shows that a sequence would be irregular if certain element is much 
larger than the others. 

Example 2.5: Consider a sequence {ci}f^^ = {10,10,1,1} and d = 2 or 3. When d = 2, we have 10^ < 1/2(10^ + 
10^ + 1 + 1). Hence {cj^^i is regular with respect to d = 2. When d = 3, we have 10^ > 1/3(10^ + 10^ + 1 + 1), 
10^ > 1/2(10^ + 1 + 1) and 1 < 1/1(1 + 1). Hence {c^}j^^ is irregular with respect to d 3 and the irregularity is fco = 2. 
This example shows that a sequence might be regular with respect to one integer but irregular with respect to another. 

The minimizers of the frame potential in ([3]) are characterized by the following lemma |19J, which will be used to prove 
the necessary and sufficient conditions for optimal placement. 

Lemma 2.6: In E"*, given a positive non-increasing sequence {ci}"^^ with uregularity as fco, if the norms of the frame 
{ipi}2=i are specified as \\(pi\\ = Ci for alH G {1, . . . , n}, any minimizer of the frame potential in ([3]l is of the form 

W^}?=l^W^}tl^W^}7=k„ + l, 

where {iPi}iZi is an orthogonal set, and {(Pij^Lfco+i is a tight frame in the orthogonal complement of the span of {^PiyiZi- 
Any local minimizer is also a global minimizer. 

a minimizer of the frame potential consists of an orthogonal set {<Pi}*5,i and a tight frame {v'iliLfco+i- 



From Lemma 



2.6 



The partition of the two sets is determined by the irregularity of the specified frame norms {ci}"^i. When the irregularity 
fco = 0, it is clear that a minimizer of the frame potential is a tight frame. As a corollary of Lemma |2.6[ the following result 
IT9) gives the existence condition of the solutions to 
Lemma 2.7: In W^, given a positive sequence {q}"^]^, there exists a tight frame {<y9i}"^i with = Ci for alH e {1, . . . , n} 
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TABLE I: Measurement models and FIMs of the three sensor types. 



Sensor type Measurement model FIM Coefficient Optimality criterion 
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solving Q if and only if {q}"^]^ is regular. 

III. Problem Formulation 

Consider one target and n sensors in M'* (rf = 2 or 3 and n > d). As shown in Table |l] we consider three types of sensors: 
bearing-only, range-only, and RSS. Suppose n sensors involve only one sensor type. Sensor networks with mixed sensor types 
are not addressed. Following ijsj, l|7)-||9], a rough estimation p e of the target position is assumed to be obtained in 
other ways. Since p is the only available information of the target, the optimal placement will be determined based on this 
estimation. As illustrated in |9|, the rough estimation of the target position can be obtained from a Kalman filter in practice and 
the optimal placement can be applied to improve target tracking performance. It should be noted that this paper only focuses 
on determining optimal sensor placements and will not discuss their application to target localization or tracking. Denote the 
position of sensor i as e M^, i G {1, . . . , n}. Then = Si — p denotes the position of sensor i relative to the target. The 
relative sensor-target placement can be fully described by {r^}"^]^. Our aim is to determine the optimal {ri}^^i such that 
certain objective function can be optimized. The distance between sensor i and the target is given by ||ri||. The unit-length 
vector gi — / 1 1 1 1 represents the orientation of sensor i relative to the target. 

A. Sensor Measurement Model and FIM 

For any sensor type in Table |l] the measurement model of sensor i is expressed as 

Zi = hi{ri) + Vi, 

where Zi E M'" denotes the measurement of sensor i, the function hi{ri) : M'" is determined by the type of the 

sensor, and Vi £ M™ is the additive measurement noise. We assume Vi to be a zero-mean Gaussian noise with covariance as 
Si = aflm G M™^™, where I„i denotes the m x m identity matrix. By further assuming the measurement noises of different 
sensors are uncorrelated, the FIM given by n sensors is expressed as 




where dhi/dp denotes the Jacobian of hi{ri) — hi{si — p) with respect to p. We refer to ||5j, ||7|-||9], pT) for a detailed 
derivation of the FIM formula in (|6|. 

The measurement models of bearing-only, range-only, and RSS sensors are given in Table |l] The measurement of a bearing- 
only sensor is conventionally modeled as one angle (azimuth) in 2D or two angles (azimuth and altitude) in 3D. The drawback 
of this kind of model is that the model complexity increases dramatically as the dimension increases. As shown in Table |lj we 
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model the measurement of a bearing-only sensor as a unit-length vector pointing from the target to the sensor A unit-length 
vector essentially characterizes a bearing and is very suitable to represent a bearing-only measurement. This model, which 
was proposed in our previous work p4) , enables us to easily formulate optimal bearing-only placement in and M^. The 
measurement model of range-only sensors given in Table [I] is the same as the one in |5|. The measurement model of RSS 
sensors given in Table |l] is a modified version of the one in ||7). Without loss of generality, we simplify the model in [Tj by 
omitting certain additive and multiplicative constants. 

The FIMs of the three sensor types are shown in Table |l] The FIM can be calculated by substituting h{ri) in to (|6]l. The 
calculation is straightforward and omitted here. As will be shown later, the coefficients {ci}"^]^ in the FIMs fully determine 
the optimal placements. Following ||5), ||7|-|[9|, we assume the coefficient to be arbitrary but fixed, (i) For bearing-only or 
RSS sensors, because q = l/(CTi||ri||), both ct^ and ||ri|| are assumed to be fixed. Otherwise, if ||rj|| is unconstrained, the 
placement will be optimal when ||ri|| approaches zero. To avoid this trivial solution, it is reasonable to assume ||ri|| to be 
fixed, (ii) For range-only sensors, because Ci = l/di, only ai is assumed to be fixed. Hence ||ri|| will have no influence on 
the optimality of the placement for range-only sensors. 

To end this subsection, we would Uke to point out that the FIMs given in Table |l] are consistent with the ones given in ||5], 
|j7j-|j9j in 2D cases. To verify that, we can substitute gi — [cos 6', , sin 6*,]^ e into the FIMs in Table |l] 

B. New Criterion for Optimal Placement 

The conventional criterion for optimal placement is to maximize the determinant of the FIM, i.e., det F. But det F is hardly 
analytically tractable in M?. In order to analytically characterize optimal placements in and M^', we will consider a new 
objective function. Let be the eigenvalues of F and A = V'^X^^Li ^i- The new objective function considered in this 

paper is \\F ~ XI dW^- Compared to the conventional one dctF, the new objective function — XI dW^ is of strong analytical 
tractability. We next formally state the optimal sensor placement problem that we are going to solve. 

Problem 3.1: Consider one target and n sensors in M'' (d = 2 or 3 and n > d). The sensors involve only one of the three 
sensor types in Table |l] Given arbitrary but fixed positive coefficients {ci}"^i, find the optimal placement such that 

{9n?=i= argmin \\F - XId\\^ (7) 

where S*^^^ denotes the unit sphere in M"^. 

Remark 3.2: The relative sensor-target placement can be fully described by {r^}"^]^. Recall ||ri|| is assumed to be fixed for 
bearing-only or RSS sensors, and ||rj|| has no effect on the placement optimality for range-only sensors. Thus for any sensor 
type, the optimal sensor placement can also be fully described by {gi}f^i. That means we only need to determine the optimal 
relative sensor-target bearings {g*}2^i to obtain the optimal placement. 

Although the FIMs of different sensor types may have different formulas as shown in Table |l] the following result shows 
that substituting the FIMs of the three sensor types into Q will lead to an identical objective function. The following result 
is important because it enables us to unify the formulations of optimal placement for the three sensor types. 

Lemma 3.3: Consider one target and n sensors in M'' (li = 2 or 3 and n > d). The sensors involve only one of the three 
sensor types in Table |l] The problem defined in Q is equivalent to 

{<?ar=i= argmin \\G\\\ (8) 
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where G = cjg^gf . 

Proof: If all sensors are bearing-only, the FIM is = ~ Sifff) ^nd then A = l/c^X^j^i-^j = ti F/d = 

(d-l)/dELicf. Hence 

z^l z^l 
n ^ n 

1=1 i=l 

If all sensors are range-only or RSS, the FIM is F = c^9i9l ™d then A — l/rfX]j=i '^j = tr F/d — l/dJ27=i ^1- 

Hence 



1=1 i=l 

Therefore, for any one of the three sensor types given in Table [I] the new objective function can be rewritten as 

n 1 ^ 




(9) 

Because l/d(X]r=i "^i)^ constant, minimizing ||G|p is equivalent to minimizing \\F — A/^lp. ■ 
One primary task of this paper is to solve the parameter optimization problem ([8]). It should be noted that we must clearly 
know the sensor type that we work with, such that the coefficients {c;}"^]^ in G can be calculated correctly according to the 
sensor type. Once {q}"^]^ have been calculated, the sensor types will be transparent to us. As a consequence, the analysis of 
optimal sensor placement in the sequel of the paper will apply to all the three sensor types. 

C. Relationship between the New and Conventional Criterions 

The new criterion for optimal placement is to minimize \\F ~ AJd|p, while the widely adopted conventional one is to 
maximize det F. We will show that the new criterion has a very close connection to the conventional one. 
Lemma 3.4: For any one of the three sensor types given in Table |l] we have 

detF < A'', 

where the equality holds if and only if 

\\F~~XIaf = 0. (10) 

Proof: For any one of the three sensor types, the FIM F is symmetric positive (semi) definite. Hence Xj is real and 
nonnegative. From the FIMs shown in Table |lj we have X]j=i -^i = trF = (d— 1) X]r=i "^1 ^^i" bearing-only sensors, and 
= tri^ = X]"=i for range-only or RSS sensors. Note {ci}^^i is assumed to be fixed. Hence X]j=i '^j is an 
invariant quantity for any one of the three sensor types. By the inequality of arithmetic and geometric means, the conventional 
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objective function dct F satisfies 

d I ^ d ^ 

" ^ (11) 



3=1 \ J=l / 



where the equality holds if and only if Aj = A for all j e {1, . . . , d}, which means 

F=~XId^\\F-~XIdf = 0. (12) 
To sum up, dcti^ is maximized to its upper bound A'' if and only if ||F — A/^IP = 0. ■ 



Based on Lemma 3.4 we next further examine the relationship between the new and conventional criterions case by case, 
(i) In R^, we have dctF = l/2((tr J')^ - tr (F^)) = 1/2 {AX^ - \\F\\^) and \\F - XhW^ = tr {F - Xh)^ = \\F\\^ - 2X^, 
which suggest 

WF-Xhf = -2detF + 2p. 

Because 2A^ is constant, minimizing — A/2IP is rigorously equivalent to maximizing det F in M^. Thus the new criterion is 
rigorously equivalent to the conventional one for all 2D cases. As a consequence, our analysis based on the new criterion will 
be consistent with the 2D results in js], l|7)-||9j. (ii) In M.^, if ||F — A/3IP is able to achieve zero, then det F can be maximized 



to its upper bound as shown in Lemma 3.4 In this case the new criterion is still rigorously equivalent to the conventional 
one. (iii) In R^, \\F ~ XI^W^ is not able to reach zero in certain irregular cases (see Section 



IV 



for the formal definition of 



irregular). In these cases detF and — A/3IP may not be optimized simultaneously. But as will be shown later, the analysis 
of irregular cases in based on the new criterion is a also reasonable extension of the analysis of irregular cases in M^. 

D. Equivalent Placements 

Before solving (|8]l, we identify a group of placements that result in the same value of ||G|p. 

Proposition 3.5: The objective function \\G\\^ is invariant to the sign of gi for all i e {l,...,n} and any orthogonal 
transformations over {gi}2^i- 

Proof: First, gigf — {—giji—gi)^ for all i E {!,..., n}, hence ||G|p is invariant to the sign of gi. Second, let 
U G R"^^"^ be an orthogonal matrix satisfying U'^U = Id- Applying U to {gi}f^i yields {g'^ = UgiYl^^. Then we 
have G' = Y.'l=ich[{9'iV = Y.l=icl{U 9i){U g^Y = UGU^ . Since G and G" are both symmetric, we have ||G'||2 = 
tr (UGU^UGU^) = tr (G^) = ||G||2. ■ 
Geometrically speaking, changing the sign of gi means flipping sensor i about the target, and an orthogonal transformation 
represents a rotation, reflection or both combined operation over all sensors. It is noticed that the invariance to the sign change 



of gi was originally recognized in ||8j for bearing-only sensors. By Proposition 3.5 we define the following equivalence 
relationship. 

Definition 3.6 (Equivalent placements): Given arbitrary but fixed coefficients {c^}"^]^, two placements {gi}"=i ^^'^ {9i}?=i 
are called equivalent if they are differed by indices permutation, flipping any sensors about the target, or any global rotation, 
reflection or both combined over all sensors. 

Due to the equivalence, there always exist an infinite number of equivalent optimal placements minimizing ||G|p. If two 
optimal placements are equivalent, they lead to the same objective function value. But the converse statement is not true in 



^3 



(b) 





Fig. 1: Examples of equivalent placements (d = 2,n = 3): (a) Original placement, (b) Rotate all sensors about the target 60 degrees clockwise, (c) Reflect 
all sensors about the vertical axis, (d) Flipping the sensor S3 about the target. 



general. In Section |V-C| we will give the condition under which the converse is true. Examples of 2D equivalent placements 
are given in Figure [T] 

IV. Necessary and Sufficient Conditions for Optimal Placement 
In this section, we will prove the necessary and sufficient conditions for optimal placements solving (|8]l. Recall G = 
Er=i '^i9^9I■ Then we have 

n n 

i=i j=i 

n 71 

1=1 j=i 

where cpi — CiQi for i e {1, . . . , n}. The vectors form a frame in W^. The objective function ||G|p is exactly the 

frame potential of the frame {^PiYl^i as shown in ([3]l. The matrix G is the frame operator Note here = c^. The coefficient 
sequence {ci\'2^i will fully determine the minimizers of UGlp. According to the irregularity of {cijf^j^, optimal placements 
are categorized as regular and irregular as shown below. 

When {ci}"^]^ is regular, the necessary and sufficient condition of optimal placement is given below. The 2D version of the 
following result has been proposed in |[5), ||7j, |[8|. 

Theorem 4.1 (Regular optimal placement): In with d = 2 or 3, if the positive coefficient sequence {c^}"^]^ is regular, 
then the objective function ||G|p satisfies 

2 



1G|P>^(Ec?) . (13) 



The lower bound of \\G\\ is achieved if and only if 



71 1 ^ 



d 



Proof: Let be the eigenvalues of G. Then Ei=i A'j = trG = ^"=1 "^f constant. Let p, = ^/d^'^^-^^ iij 
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l/'^X]"=i cf- It is obvious that 

The lower bound of UGlP is achieved if and only if /ij — fi for all j E {!,. . . ,d}, which implies G = flld, i.e., the equation 
( [l4| i. By denoting (yCi = Cigi, ( [T4] i becomes '/'ivf = "^f^d which is the same as (|2]i. Thus a regular optimal 

placement solving ^14) actually corresponds to a tight frame. Because {cj}"^j^ is regular, by Lemma [zT] there exist optimal 
placements solving ( [l4j l. ■ 
We call a placement regular when its coefficient sequence is regular, and regular optimal when it solves ([14]). To obtain a 
regular optimal placement, we still need to solve ([T4|. Details of the solutions to ([T4| will be given in Section V-A 



3.4 



the 



Remark 4.2: When ||G|p reaches its lower bound l/rf(X]r=i '^?)^' '^^ '^^^^ 11-^ ^ ^-^<i|P = by (|9]). By Lemma 
conventional objective function det F will be maximized to its upper bound. Therefore, regular optimal placements not only 
minimize the new objective functions ||G|p and ||F — A/(j|P but also maximize the conventional one det F in and M.^. 

When {ci}"^]^ is irregular, ( [T4] l will have no solution. Then the the necessary and sufficient condition of optimal placement 
is given below. The 2D version of the following result has been proposed in |j5), |[7], |[8). 

Theorem 4.3 (Irregular optimal placement): In M"* with d = 2 or 3, if the positive coefficient sequence {c^}"^]^ is irregular 
with irregularity as fco > 1, without loss of generality {ci}f^i can be assumed to be a non-increasing sequence, and then the 
objective function ||G|p satisfies 

fco / n \2 

ii^sE'J + j^ E - ('5) 

i=l " \'i = feo + l / 

The lower bound of ||G|p is achieved if and only if 

{9^}7=l={9^}tl^{9^}7=ko + l, (16) 

where {gi}'^li is an orthogonal set, and {ffi}iLfco+i forms a regular optimal placement in the (d — fco) -dimensional orthogonal 
complement of {gi}'lZi- 



Proof: Recall ||G|p is the frame potential of the frame {<y5i}"=i where ipi — Cigi. From Lemma 2.6 the minimizer of ||G|p 
is of the following form: {cigi}^!^^ is an orthogonal set, and {cigi}'^^i^^^i is a tight frame (i.e., a regular optimal placement) 
in the orthogonal complement of {ci^ij^^j^. 

Let $1 = € E'*><'=°, $2 = yko+i,---,Vn] e Rdx{n-ko)^ $ ^ [$i,$2] G M'^><". When {5,}^^, is of the 

form in ([T6]l, the columns of $1 are orthogonal to those of $2- Then 

||G|p = tr = tr (^f + tr ($^«>2)^ 

Because {gi}'^li is an orthogonal set, we have tr (<i>|"<i>i)^ = X^i^i llv'ill^ = X]i^=i '^f ■ Because {gi}7^f.^-^^i is a regular optimal 



placement in a (d — fco) -dimensional subspace, we have tr {^2^2)^ = ^/{d — ^o){J2"=ko+i "^i )^ Theorem 4.1 Therefore, 
when {gi}f^i is of the form in ([T6|, the objective function ||G|p reaches its lower bound as shown in ( [T5| l. ■ 
We call a placement irregular when its coefficient sequence is irregular, and irregular optimal when it is of the form in 
because {gi}2=ko+i is a regular optimal placement in a [d— fco) -dimensional space, an irregular optimal 



rSll. In Theorem 



4.3 



placement problem will eventually convert to a regular^ one in a lower dimensional subspace. 
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(c) 

Fig. 2: An illustration of the three kinds of irregular optimal placements in and R"^. (a) d = 2, kg = I; (b) d = 3, fco = 1; (c) d = 3, fco = 2. 

Apparently the irregularity of {c^}"^]^ plays a key role in determining optimal placements. Recall the irregularity fcp of an 
irregular sequence with respect to d satisfies 1 < fco < d — 1. In M^, we have d = 2 and hence fco = 1; in we have d = 3 
and hence fco = 1 or 2. Thus there exist only three kinds of irregular optimal placements in and R^. By Theorem 4.3 
these three kinds of irregular optimal placements can be intuitively described as below. 

(i) Irregular optimal placement in with irregularity fcg = 1- the vector gi is orthogonal to {.gi}"^2' ^"'^ {5i}"=2 
coUinear. See an illustration in Figure |2] (a). 

(ii) Irregular optimal placement in M.^ with irregularity fcp = 1: the vector gi is orthogonal to {5i}r=2' {9i}?=2 form a 
regular optimal placement in the 2D plane perpendicular to gi. See an illustration in Figure |2] (b). 

(iii) Irregular optimal placement in M.^ with irregularity fco — 2: the vectors gi, 52 and {gi}f^^ are mutually orthogonal, and 
{9i}?=3 are collinear See an illustration in Figure 1(c). 

We now intuitively explain why irregular optimal placements are the way as described above. The sequence {cijf^j^ is 
irregular when certain Ci is much larger than the others. The coefficient Ci actually is the weight for sensor i. The larger the 
weight Ci is, the more sensor i contributes to the FIM. Recall q = l/((Ti||ri||) or q = l/(Ji- If sensor i is very close to the 
target (i.e., ||ri|| is small) and its measurement is very accurate (i.e., ai is small), then Ci will be very large and sensor i may 
dominate all the others. Thus the dominant sensor can measure the target sufficiently accurately in one dimension, and all the 
other sensors should measure the target in the orthogonal complement in order to improve the overall measurement accuracy. 

To make our analysis more general, we do not assume <t/s to be identical. It is also meaningful to check the special case, 
(Ti = aj for all i ^ j, which often arises in practice. For bearing-only or RSS-based sensors, the coefficient is c,; = l/((Ti||ri||). 
So when <Ti = aj for all i 7^ j, a regular sequence {ci}f^i implies the following equation according to the fundamental 
inequality (|5]): 

1 1 " 1 

which geometrically means no sensor is much closer to the target than the others. Conversely, if certain sensors are much 
closer to the target than the others, then the sequence {cijf^i will be irregular. The 2D version of inequahty ^TT) has been 
proposed in ||5j, |j7|, ||8|. For range-only sensors, the coefficient is Ci ~ l/oi. If at = aj for all i ^ j, then Ci = Cj and {ci}f^i 
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is regular with respect to any d < n. 

We next consider an important special case n = d, i.e., the sensor number equals to the dimension of the space. This case is 
important because the optimal placement will be independent to the coefficients {ci}"^]^ in this case. The optimal placement 
in the case of ?i = d = 2 has been solved by ||5), ||7j-||9j. 

Theorem 4.4: In M'' with d = 2 or 3, if n = d, the objective function \\G\^ satisfies 

i=l 

The lower bound of UGlP is achieved if and only if {gi}f^i is an orthogonal basis of M'*. 
Proof: Since G — X^iLi ^1.9i9i ^^'^ 9i 9i — 1 foi" * ^ {1j ■ • ■ i "■}' we have 

\\Gf^tv{G^) 

i=i j=i 

i=l j = i=l 
d 

where the equality holds if and only if gj gj = for all i.j E {1, . . . ,d} and i ^ j. ■ 



Theorem 4.4 can also be proved as a corollary of Theorem |4. 1 1 and Theorem 4.3 But as shown above, we can also directly 
prove it without using frame theory. This can be explained from the redundancy point of view. Recall the constant n/d reflects 
the redundancy of the system. There will be no redundancy when n/d = 1. Then frames are no longer necessary for the 
optimality analysis. 

V. Analytical Properties of Optimal Placements 

In this section, we further explore a number of analytical properties of optimal placements in 2D and 3D. Theorem |4.3| 
implies that an irregular optimal placement problem can be eventually converted to a regular one in a lower dimensional space. 
Hence we will only focus on regular optimal placements. 



A. Explicit Construction 

The work in ||5), Q-Q has proposed explicit construction methods for some special 2D optimal placements. The construction 
of generic optimal placements in 2D or 3D is still an open problem. From Theorem |4. 1| we know a regular optimal placement 
is essentially a tight frame. In fact, construction of tight frames in arbitrary dimensions has been extensively studied. Thus by 
referring to pO) , p2)-p4), to name a few, we are able to construct optimal placements with an arbitrary number of sensors 
and arbitrary but fixed coefficients in 2D and 3D. 

We here present a proof of the necessary and sufficient condition for 2D regular optimal placements without using frame 
theory. In the meantime, we propose an explicit construction method for arbitrary 2D regular optimal placements. The necessary 
and sufficient condition for 2D optimal placements has already been given in |j5), ||7j, |[8j, where the sufficiency, however, is 
not proved. We will prove the sufficiency by construction. The following lemma can be found in Q, ||7)-||9], |[T8|, p5| , p6) . 
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Lemma 5.1: In M^, the unit-length vector gi can be written as gi = [cos6i,smdi]'^ . Then ( [T4) l is equivalent to 



where gi = [cos 20^, sin 26*^] . 

Proof: Substituting gi — [cos 0^, sin ( 



E 

i=l 

into ( [T4] i gives 

n 



0, 



(18) 



^ cos 20^ 



29, 



i sin 26*, 

— i COS 2di 



which is equivalent to ( [T8] l. 



By Lemma 5.1 the matrix equation ([14]) is simplified to a vector equation ([T8|. In order to construct {gi}f^i solving ([14]), we 
can first construct {gi}2=i solving ( [TS] ). Once g,j = [cos 20^, sin 26*^]-^ is obtained, gi can be retrieved as gi = ±[cos sin6'i]-^. 



Note the sign changes of gi give equivalent optimal placement as mentioned in Definition 3.6 

Theorem 5.2: In M^, given a positive sequence {q}"^]^, there exists {gi]^^i with \\gi\\ = 1 solving ([TSjl if and only if 



max c < 

j = l,...,Tl 



1 " 

E 



cl 



(19) 



E 



Proof: Necessity: If J2'i=i ^Idi = 0, then c^gj = J2i^j cfgi for all j e {1, ... , n}. Hence c| ||c|.gj || 
W^fdiW = Ei^^j c|. Then adding on both sides of the inequality gives 2c| < J2i=i c? ■ 
Sufficiency: If cj < l/2J27=i J ^ {1; ■ ■ ■ j*^}' it is obvious that there always exists an index no (2 < no < n) 

such that 



rio-l 



1^ 

i<2E' 

1=1 

1 " 



(20) 
(21) 



When no < n, denote 



^2 - 



(22) 



Obviously £i + £2 + ^3 = Ef^i cf ■ From ([19]), c,,^ < 1 /2 X;r=i cf and hence £1 + £3 > ^2- From ([20|, £1 < 1/2 X;r=i cf and 
hence £2 + ^3 > ^i- From ([21}, £1 + £2 > 1/2 ELi and hence ^1 + ^2 > 4- Therefore, £1, ^2 and £3 satisfy the tiiangle 
inequality and can form a triangle. Choose gi — ■ ■ ■ — (jno-i- Then ^"=7^ '^i9i — ^iffi- Choose 5^0+1 = ■ ■ ■ ~ (jn- Then 
I]"=no+i 5i = ^33n- Then ([T8| becomes 



^l5l + ^25no + 45n = 0. 



(23) 



We can choose gi,gno and that align with the three sides of the triangle with side length as £1, £2 and £3, respectively (see 
Figure [3}. Then ( [23] ) and consequently ( [TS] ) can be solved. When no = n, the above proof is still valid. In this case, we have 



and £1 — £2, and ([23]) becomes gi + g. 



0. 



From the proof of Theorem 5.2 an explicit construction of 2D regular optimal placements can be summarized in Algorithm 
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gi {i > no) 


(2/ "\J3 






/ 

9no 


gi (i < no) 



Fig. 3: Geometric illustration for the 2D construction. 



Algorithm 1 Construction of 2D regular optimal placements {gi}f^i with coefficients {ci}"^]^. 
1: Choose no satisfying ( |20] i and pT| ). Then compute ii, £2 and £3 in ( |22] i. 

2: Compute interior angles 0112 and of the triangle with side lengths as £1, £2 and £3 (See Figure [3|. 
3: Choose gi = [1,0]'^ for i £ {l,...,no - 1}, g„o = [cos((7r + ai2)/2), sin((7r + Q!i2)/2)]^, and g^ = 
[cos((7r - ai3)/2), sin((7r - ai3)/2)]^ for i G {no + 1, . . . , n}. 



[T] The following example illustrates the construction method in Algorithm [T] 

Example 5.3: In M^, given six bearing-only sensors with sensor-target ranges respectively as ||ri|| ~ 5, ||r2|| ~ 6, ||r3|| = 7, 
Wr^W = 8, llrsll = 9, and \\rg\\ = 10. The measurement noise variance is cr.j ~ 1 for all i e {1, . . . ,6}. Recall Ci = l/{(Ti\\ri\\) 
for bearing-only sensors. Then cj = 0.0400, 4 = 0.0278, cj = 0.0204, = 0.0156, = 0.0123, cj = 0.0100, and 
l/2^^^j^cf = 0.0631. It is easy to check the sequence {q}^^]^ is regular. Because < l/2X]i=icf and c\ + > 
l/2ELicf, choose no = 2. Hence £1 = 0.0400, £2 = 0.0278, and £3 = 0.0584. Then ai2 = 2.0560rad and ais = 
0.4344 rad. As instructed in Algorithm [l] choose gi = [1,0]^, 32 = [0.8563, -0.5165]'^, 33 = ■ ■ • = 56 = [0.2155,0.9765]^. 
It can be verified that J2Li ^Ig^al = 1/2 ELi cf ^2. 



B. Equally-weighted Optimal Placements 

The coefficient Ci actually is the weight of sensor i. A placement is called equally-weighted if ci = • • • = €„. Since {ci}^^i 
is regular with respect to any d < n if ci • • • = c„, equally-weighted placements must be regular. For bearing-only or RSS 
sensors, equally-weighted means di — aj and ||r,H = ||rj|| for all i 7^ j as — l/{<Ti\\ri\\). The corresponding geometry is 
that all sensors are restricted on a 2D circle or a 3D sphere centered at the target. For range-only sensors, equally-weighted 
means ct; = <Tj for all i ^ j as Ci — l/(Ji. 

Equally-weighted placements are important because they often arise in practice and have some important special properties. In 
the equally-weighted case, ( [T4] l is simplified to ^"=1 9i9i ~ n/did, which implies that an equally-weighted optimal placement 
is essentially a unit-norm tight frame fls) , | |26) . In M^, an equally-weighted placement is optimal if n (n > 3) sensors are 
located at the vertices of an n-side regular polygon ijsj, 0-|j9), 1 18 1, |j26) as shown in Figure |4] In M'^, an equally-weighted 



placement is optimal if n sensors are located at the vertices of a Platonic solid 1 18 1, |26|. There are only five Platonic solids as 
shown in Figure |5] It should be noted that equally-weighted optimal placements are not limited to regular polygons or Platonic 
solids. In Section |V-D we will show more examples of equally-weighted optimal placements. 
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(a) (b) (c) (d) (e) 

Fig. 4: Examples of 2D equally-weighted optimal placements: regular polygons. Red square: target; blue dots: sensors. 




(a) (b) (c) (d) (e) 

Fig. 5: Examples of 3D equally-weighted optimal placements: Platonic solids. Red square: target; blue dots: sensors, (a) Tetrahedron, n = 4. (b) Octahedron, 
n = 6. (c) Hexahedron, n = 8. (d) Icosahedron, n = 12. (e) Dodecahedron, n = 20. 



C. Uniqueness 

Due to placement equivalence, there exist at least an infinite number of equivalent optimal placements minimizing |lG|p. 
It is interesting to ask whether all optimal placements are mutually equivalent. Now we give the condition under which all 
optimal placements are mutually equivalent, or in other words, the optimal placement is unique up to the equivalence. Again 
we only consider regular cases. 



According to Theorem 4.4 it is clear that the optimal placement is unique in the case of n = d. We next show the regular 
optimal placement is also unique in the case of n = d + 1 (i.e., three sensors in M? or four sensors in M.^). The uniqueness 
will be proved by construction, which is inspired by the work in p5] on unit-norm tight frames. 

Theorem 5.4: In M'* with d = 2 or 3, if n = d + 1, given a regular coefficient sequence {ci}fl^^, the regular optimal 
placement {gi}i^i is unique up to the equivalence in Definition 



3.6 



Proof: Suppose {gi}fl^i is a regular optimal placement solving ( [l4| i. Denote ipi = Cigi and $ = [tpi, . . . ,tpd+i\ & 
|dx(d+i) jjjgjj (jY4j| (.^jj written in matrix form as <I><I>-^ = l/c^X^itti cf/d- Hence $ has mutually orthogonal rows with 



row norm as y l/rfX]f=i "^f- Let x = [xi,. .. ,Xd+i\ G M''+^ be a vector in the orthogonal complement of the row space 
of $. Assume ||a:;|| = \J ^1- Adding x^ after the last row of $ yields an augmented matrix ^'aug = ['^■"^,2^]^ G 

]]j(c!+i)x(ci+i) jj. clear that $aug"J'rug — 5 Si^i '^f-^d+i- Thus $aug is a scaled orthogonal matrix and its columns are 
mutually orthogonal. The jth column of $aug is VpJ i Xj]^ E M''+^ for all j G {1, . . . , rf + 1}. Note the column norm of <&aug 
^l/dJ2i=i cf ■ Then we have WipjW^ + x'j — l/dJ2i=i c| and hence 



± 



, d+l 



i=l 



The regularity condition ensures l/^X^i^i "^1 ~ "^f — ^ 
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Algorithm 2 Construction of the unique regular optimal placement {gijfjli with coefficients {ci} 



id+l 



1: Choose X = [xi, . . .,Xd+i] E with xj = i^l/d^^tti cf - cj for i E {1, . . . ,d + I}. 

2: Use the singular value decomposition (SVD) to numerically compute an orthogonal basis of the orthogonal complement 
of X. Let X = UY.V'^ be an SVD of x, where U G ]R(rf+i)x(ci+i) is an orthogonal matrix. 



(27) 



3: Let Ui denote the ith column of U. Then x = "^f^i' ^'^'^ ^ '-^'^ constructed as 



^ d+l 



Ddx(d+1) 



4: Compute gi = (pi/ci for i <E {1, . . . , + 1}. 



By reversing the above proof, we can obtain an explicit construction algorithm for optimal placement with n — d + 1 as 
shown in Algorithm |2] The rest is to prove the constructed optimal placements are mutually equivalent. First, given a vector 
satisfying (|24]|, let $ and $' be two different bases of the orthogonal complement of x. Due to orthogonality, there 



X e 



exists an orthogonal matrix U G such that 

U 

Write U as 

U = 









T 






X 




X^ 



(25) 



U21 U22 



(26) 



where Un £ R'^'"^,Ui2 € ]R''^\t/2i € R^""^, and U22 G K. Substituting ^ into (|25]) gives C/2i$ + (J722 - l)a;^ = 0. 
Since the rows of $ and x^ are linearly independent, we have U21 ~ 0, U22 = 1- Thus U12 = and C/ii$ = <!>'. Therefore, 



the placements described by $ and $' are differed only by an orthogonal transformation Un. From Definition 3.6 the two 
placements are equivalent. Second, let E e be a diagonal matrix with diagonal entries as 1 or —1. Given arbitrary 

X and x' both satisfying ( |24] l, there exists an E such that x' = Ex. Note E is also an orthogonal matrix. It can be analogously 
proved that the optimal placements would be differed by an orthogonal transformation and a number of flipping of sensors 



about the target. From Definition 3.6 these placements are also equivalent. ■ 
From the proof of Theorem |5.4[ an explicit construction of the unique regular optimal placement in the case of n = d + 1 

can be summarized in Algorithm |2] The following example illustrates the construction method in Algorithm [2] 

Example 5.5: In M.^, given four bearing-only sensors with sensor-target ranges respectively as ||ri|| = 20, ||r2|| = 21, 

Wr^W = 22 and ||r4|| = 23. The measurement noise variance of the ith sensor is ai = 0.01 with i E {1,...,4}. Recall 



l/{ai\\n\\) for bearing-only sensors. Then = 25.00, 4 = 22.68, c§ = 20.66, c| = 18.90 and l/3X;i=r 



29.08. 



The sequence {ci}j^i is regular. From ( |24] i, choose x — [2.02,2.53,2.90,3.19]^. Compute the SVD of x and use ( pTj l to 



compute $ as 



$ = 



-2.5307 4.5286 -0.9906 -1.0891 
-2.9016 -0.9906 4.2568 -1.2487 
-3.1901 -1.0891 -1.2487 4.0197 



It can be verified ^ti cfftsf = <^<^^ = 1/3 Eti ^fh- 
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(a) (b) 

Fig. 6: The unique equally-weighted optimal placements with n = 3 in R^. Red square: target; blue dots: sensors, (a) Regular triangle, (b) Flip si about the 
target. 




(a) 

Fig. 7: The unique equally-weighted optimal placements with n = 4 in ] 
the target, (c) Flip S4 and S3 about the target. 



(b) 



(c) 



Red square: target; blue dots: sensors, (a) Regular tetrahedron, (b) Flip S4 about 



Figure [6] and Figure [7] show examples of unique optimal placements. From the last subsection, we know a regular triangle 



and a regular tetrahedron are equally-weighted optimal. By Theorem 5.4 they are also unique up to equivalence. Hence the 
two equivalent placements in Figure [6] consist of all possible forms of the equally-weight optimal placements with n = 3 in 
M^. The three equivalent placements in Figure [t] consist of all possible forms of the equally-weight optimal placements with 

n = 4 in M.^. 

When n > d+1, the regular optimal placement may not be unique. In the next subsection, we will give examples to show 
the optimal placement may not be unique when n > 4 in or n > 6 in M"^. Now a question remains: whether the regular 
optimal placement with n = 5 in M.^ is unique up to the equivalence. The answer is negative. The following is an explanation 
and an explicit construction of the regular optimal placement with ?i = 5 in M"^. 



Suppose the sequence {cijf^j^ is regular with respect to = 3. Denote (pi = Cigi and $ = [(^1, . . . , ip^] G 



p3x5 



Then ^ 



becomes <i><i>-'" — 1/3 c|-?^3. There always exists $' 
space of $ such that 



p2x5 



in the orthogonal complement of the row 



IT 



1^ 



= l/3Ei=ic| and = l/3Ei=ic|/2- Thus represents a 2D regular optimal 



which implies HiyajlP 

placement with = \J^/iYTi=i '^f ^ '^^ all j G {1, ■ • ■ ,5} (it can be verified {||<y3j||}^^]^ is regular with respect to 

d = 2). Therefore, to obtain $, we can first construct {<Pj}j^i using Algorithm [T| for example, and then find $ in the orthogonal 
complement of the row space of Since {p'i}^^i may have non-equivalent solutions, {</3i}f^i would not be unique up to 
the equivalence. 
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D. Distributed Construction 

When there are a large number of sensors, it might be inconvenient to design the optimal placement involving all sensors. 
The following property can be applied to construct large-scale optimal placements in a distributed manner. The 2D versions 
of the following result have been proposed in Q, |]7], |[8|. 

Theorem 5.6: The union of multiple disjoint regular optimal placements in {d ~ 2 or 3) is still a regular optimal 
placement in M'^. 

Proof: In W^, consider multiple disjoint regular optimal placements: {ci,gi\i^i^ with Ik as the index set of the fcth 
placement (fc = 1, . . . ,q). The term disjoint as used here means that different placements share no common sensors. Define 
I = {Jl^i^k- If I ■ I denotes the cardinality of a set, then \I\ — J2k=i 



For the fcth placement, since {ci, gi}i^x^ is regular optimal in R'^, from Theorem 4.1 we have 



For the union placement, we have 



9 

j€l k=l ielfc 



k—1 i^Xk 



By Theorem 4.1 the union placement is regular optimal in M . ■ 
Theorem 5.6 implies that a large-scale regular optimal placement can be constructed in a distributed manner: firstly divide 
the large-scale placement into a number of disjoint regular sub-placements, secondly construct each regular optimal sub- 
placement, and finally combine these optimal sub-placements together to obtain a large regular optimal placement. We call this 
kind of method distributed construction. Because the combination of the optimal sub-placements can be arbitrary, distributed 
construction will lead to an infinite number of optimal placements for the large system. These optimal placements have the 
same FIM and ||G|p, but they are generally non-equivalent. From Theorem 5.6 we also know only regular placements can 
be possibly divided into some regular subsets. 

The distributed construction method is very suitable for constructing equally-weighted optimal placements. That is because 
any equally-weighted sequence {ci}^!^^ is regular if fc > d, and thus {c;}"^]^ can be easily divided into a number of regular 
subsets. We next present two instances, (i) For any integer n > 4, it is obvious that there exist nonnegative integers mi and 
777,2 such that 77 can be decomposed as ri = 2?77i +87772. Thus in M? we can always distributedly construct an equally- weighted 
optimal placement with 77 > 4 by using the ones with 77 = 2 or 3. (ii) For any integer 77 > 6, there exist nonnegative integers 
7771, m2 and 7773 such that n = 3mi + 47772 + 5777,3. Thus in we can always distributedly construct an equally- weighted 
optimal placement with 77 > 6 by using the ones with 77 = 3, 4 or 5. As a consequence, noticing distributed construction 
yields an infinite number of non-equivalent optimal placements, equally-weighted placements with 77 > 4 in M'^ or 77 > 6 in 
M'^ always have an infinite number of non-equivalent optimal solutions. 

We show some typical examples in Figure [8] to illustrate the distributed construction. From Section V-B[ we know a regular 
triangle or regular tetrahedron placement is equally-weighted optimal. The placement shown in Figure [8] (a), (b) or (c) is 
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(d) (e) 
Fig. 8: Examples of optimal placement by distributed construction. Red square: target; dots: sensors. 



distributedly constructed by two regular optimal placements with n = 3 (the sensors with the same color form a triangle 
optimal placement). The placement shown in Figure [8] (d) or (e) is distributedly constructed by the optimal placement with 
n = 4 as shown in Figure |7] (c), which is equivalent to the regular tetrahedron. Thus by Theorem |5.6| all placements in Figure 
|8] are regular optimal. 

VI. Numerical Verification 

In order to verify our previous analytical analysis, in this section we solve the parameter optimization problem ([8| from a 
numerical perspective and then present some numerical simulations. More specifically, we employ Lyapunov approaches to 
design a centralized gradient control law which can numerically minimize the objective function HGH^ given an appropriate 
initial point. The control law can be applied to numerically construct generic regular and irregular optimal placements in 2D 
and 3D. 

Assume the motion model of sensor i to be Si = Ui, where Ui e is the control input. Then we have fi = Ui because 
ri = Si — p and the target position estimation p is given. Let r = [rf, . . . , r^Y' ^ M''". Denote /3 as the constant lower bound 
of ||G|p as shown in Theorem 
the Lyapunov function as V{r) ~ l/4(||G|p — Clearly V is positive definite with respect to £q. Denote dV/dri as the 
Jacobian of V with respect to r^. Then we have 

" r>V " r2 

where Pi — Id ~ gigj is an orthogonal projection matrix satisfying — Pi, Pf — Pi, and Null(Pi) — span{(7i}. Null(-) 
denotes the nuU space of a matrix. Design the gradient control law as 



4.1 



and 



4.3 Then the optimal placement set is £q = {r E 



|G|p - /? = 0}. Choose 



ri = -PiGgi, 



(28) 
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Fig. 9: Gradient control of equally-weighted (regular) placements with n = 4 in ' 
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(a) n = 3, ko = 1 (h) n = A, ko = I (c) n = 4, ko = 2 

Fig. 10: Gradient control of irregular placements in M"^. 



such that 

n 2 

and V = when PiGgi = for alH e {1, . . . , n}. 

Proposition 6.1: For any initial condition r(0) G M*' with ||r,;(0)|| 7^ for all i E {1, . . . ,n}, the solution to the nonlinear 
r-dynamics ( |28] l asymptotically converges to the set 

£ = {r eS : P,Gg^^O, i = l,...,n}, 

where S = {r e M*' : \\n\\ = ||r,(0)||, i = l,...,n}. 
Proof: The time derivative of ||ri|| is 

^ = {^h = ~9lP^Gg, = 0. (29) 
The last equality uses the fact gfPi = 0. By ( |29| ) we have = ||ri(0)|l ^ 0. Hence 5 is a positive invariant set with 
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respect to the r-dynamics. The set S consists of a group of spheres in M.'^ centered at the origin. Thus S is compact. Note 
V ^ and f = for all points in £. By the invariance principle ( jz?] , every solution starting in S asymptotically converges 
to £. m 



By Proposition 6.1 the r-dynamics converge either to the optimal placement So or the set £ \ Eq- By introducing Lagrange 
multipliers ji, i = 1,. .. ,n, the constrained optimization problem ^ is equivalent to minimizing the Lagrangian function 
L ~ ||G|p + Y^^=ilii9l 9i ~ !)■ From dL/dgi = 0, a short calculation shows that £ is the critical point set, which 
consists of not only minimizers of ||G|p (i.e., optimal placements) but also saddle points and maximizers of ||G||^ (i.e., non- 
optimal placements). The sets £o ™d £ are equilibrium manifolds. It is noticed that nonlinear stabilization problems involving 



equilibrium manifolds also emerge in formation control area recently |28|-|30|. It is possible to conduct strict stability analysis 
including identifying the attractive region of £o by using Center manifold theory p8| , (30| or differential geometry p9) . But 
that will be non-trivial because the geometric structure of £0 is extremely complicated as shown in | |3T| . 

Figure |9] and Figure 10 show a number of optimal placements obtained by the proposed gradient control law. Due to space 



limitations, we only show examples in 3D. The three final converged placements in Figure |9] are actually the three regular 
optimal ones shown in Figure [7] The three final placements in Figure 10 are the two as illustrated in Figure [2] (b) and (c). 



Clearly the numerical results are consistent with our previous analytical analysis. The optimality error refers to the difference 
between ||G||^ and its lower bound given in ( [T3| ) or ([TSj. The optimality error can be used as a numerical indicator to evaluate 
the optimality of a placement. As shown in Figure |9] and Figure 10 the optimality errors all converge to zero. 



A. Simulation 

More simulations will be presented to demonstrate the performance of the gradient control. Our previous analytical charac- 
terization of optimal sensor placements will be verified by the simulations. 

1) Scenario 1 for bearing-only or RSS-based sensors: For bearing-only or RSS-based sensors, the sensor-target ranges are 
assumed to be fixed. The control law (|28]l naturally fulfills this assumption. Examples of 2D and 3D optimal placements 



achieved by the gradient control are shown in Figure 11 and Figure 12 respectively. The noise covariances of all sensors are 
chosen to be identical in the simulation. 



Figure 1 1 shows 2D examples. In (a) or (e), the final angle subtended at the target by the two sensors is 90 degrees. In (b) 
and (c), the two optimal placements are equivalent. In (d), the optimal placement is a distributed construction by two optimal 
two-sensor sub-placements. The optimal placements in (f), (g), and (h) have non-equal sensor-target ranges. The placement in 
(g) is UTegular. The sensor with shortest sensor-target range drives the other two sensors to be collinear with the target. 



Figure 12 shows 3D examples. In (a) or (e), the final angle subtended at the target by any two sensors is 90 degrees. The 
final placement in (b) forms a tetrahedron. The final placements in (c) and (d) are both equivalent to the one in (b). An 
optimal placement with irregularity as fco = 1 is shown in (f). The sensor with shortest sensor-target range drives the other 



three sensors to an orthogonal plane, in which the three sensors form a 2D optimal placement like the one in Figure 1 1 (c). 
An optimal placement with irregularity as fco = 2 is shown in (g). These sensors are mutually orthogonal except the two with 
long sensor- target ranges collinear. A complicated optimal placement with ten sensors is given in (h). 

These simulation results are consistent with our previous optimality analysis. It is also observed that sensors can converge 
to an optimal placement which is close to their initial positions. 

2) Scenario 2 for range-only sensors: For range-only sensors, the sensor-target ranges are not required to be fixed. We 
consider a 3D scenario which might be applied to cooperative air and ground surveillance p2). Consider two unmanned aerial 
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S 1. 




(e) n = 3, fco = 1 (f) n = 4, fco = 1 (g) n = A, ko = 2 (h) n = 10, fco = 

Fig. 12: Examples of 3D optimal placements achieved by the gradient control law. /cq: in'egularity; green diamonds: initial sensor positions; red dots: final 
sensor positions; blue curves: sensor trajectories; black square: target. 
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(a) 3D view (b) Top view 

Fig. 13: Cooperative air and ground surveillance. Green diamonds: initial vehicle positions; red dots: final vehicle positions; blue curves: vehicle trajectories; 
black square: ground target. 
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Fig. 14: The error between ||G|P and its lower bound. 

vehicles (UAVs) and two unmanned ground vehicles (UGVs). Each vehicle carries a range-only sensor, which can measure 
the distance from the ground target to the vehicle. UAVs fly at a fixed altitude and UGVs move on the ground. In addition to 
the objective function ||G|p, we introduce the following external potential for sensor i: 

where 63 — [0, 0, 1]^ and £i is the required altitude of the vehicles. Here the external potential is used to control the sensor- 
target range to fulfill level-plane position constraints. Denote Vj — ||G||^ as the inter-sensor potential, and Ve = ^siri) 
as the external potential of all sensors. The total potential function is V ^Vj + Ve- And the gradient control law becomes 
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where V^^ is the gradient with respect to r^. The external potential Ve should be compatible with the inter-sensor potential Vj, 
otherwise, the external potential will affect the convergence to optimal placements. By compatible, we mean Vf+V^j = 
if only if V^jV/ — and V,-. Ve = 0. Since Vr^Vj _L r.^, Ve is compatible with Vj if there is no make Vr^VE -L ^i. One 
compatible example is that all range-only sensors move on the boundary of a convex area containing the target |[9). 

In the simulation, we choose ^1=^2 = 10 and £3 — £4 = 0. The noise variances of different sensors are the same. Figure 



13 shows the trajectories and the converged placement of the vehicles. For compUcated placements, it is usually difficult to 
intuitively judge their optimality. We can evaluate the optimality by using the error between UGlP and its lower bound given 
in ( [T5| ) and ^T3\ . As shown in Figure [14] the optimality error converges to zero. 

VII. Conclusions 

The main contribution of this paper is that we present a unified way to analytically characterize optimal placements of 
bearing-only, range-only, or RSS sensors in 2D and 3D. The necessary and sufficient conditions for optimal placement in 2D 
and 3D are proved. The gradient control proposed in this paper is already sufficient for the purpose of numerical verification 
of the analytical analysis. But it is a centralized control based on all-to-all communications and ensures local stability of the 
optimal placement set. It is meaningful to study distributed or global stability guaranteed control laws or numerical algorithms 
in the future. 
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